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Abstract 

The phase diagram of strongly interacting matter at non-zero temperature and 
baryon chemical potential is calculated within a 3-flavor NJL-type quark model with 
realistic quark masses. The model exhibits spontaneous chiral symmetry breaking 
as well as diquark condensation in the two-flavor color-superconducting phase and 
in the color-flavor locked phase. We investigate the color-flavor unlocking phase 
transition, taking into account self-consistently calculated effective quark masses. 
We find that it is mainly triggered by a first order phase transition with respect 
to the strange quark mass. It takes place at much higher values of the chemical 
potential than the transition to the hadronic phase such that we find a relatively 
large region in the phase diagram where the two-flavor color-superconductor seems 
to be the most favored state. 



1 Introduction 



The structure of the QCD phase diagram is one of the most exciting topics in the field 
of strong interactions (For review see, e.g. ]T|, [| [|). For a long time the discussion 
was restricted to two phases: the hadronic phase and the quark-gluon plasma (QGP). 
The former contains "our" world, where quarks and gluons are confined to color-neutral 
hadrons and chiral symmetry is spontaneously broken due to the presence of a non- 
vanishing quark condensate <ft — (V'V'}- m t ne QGP quarks and gluons are deconfined and 
chiral symmetry is (almost) restored, 0^0. 

Although color-superconducting phases were discussed already in the '70s [§, |6j and 
'80s J7|, until quite recently not much attention was payed to this possibility. This changed 
dramatically after it was discovered that due to non-perturbative effects, the gaps which 
are related to these phases could be of the order of A ~ 100 MeV |8|, much larger than 
expected from the early perturbative estimates. Since in standard weak-coupling BCS 



theory the critical temperature is given by T c ~ 0.57 A(T = 0) fllPl , this also implies a 



sizable extension of the color-superconducting phases into the temperature direction fil l. 



Hence, color-superconducting phases could be relevant for neutron stars [12], |13[ and - in 
very optimistic cases - even for heavy-ion collisions JR]. 

Many different diquark condensates are allowed by Pauli principle 0. However, most 
interactions, like one-gluon exchange or instanton mediated interactions favor a conden- 
sation in the scalar color-antitriplet channel. In general, these condensates read 

saa' = Cj 5 t a \a>iI>) , (1) 

where both, ta and \a> are the antisymmetric generators of SU(Nf) and SU(N C ), acting 
in flavor space and in color space, respectively. Throughout this paper, we will restrict 
ourselves to the physical number of colors, N c = 3. Then the \a> denote the three 
antisymmetric Gell-Mann matrices, A2, A5 and A7. 

For two flavors (Nf = 2), the flavor index in Eq. (P is restricted to A = 2, describing 
the pairing of an up quark with a down quark. The three condensates S2A', A' = 2, 5, 7, 
form a vector in color space, which always can be rotated into the A' = 2-direction. Hence 
the two-flavor superconducting state (2SC) state can be characterized by 

s 22 + and saa' = if (A, A') + (2, 2) . (2) 

Due to the presence of S22, color SU(3) is spontaneously broken down to SU(2). For 
massless up and down quarks the 2SC state is invariant under chiral SU{2)l x SU(2)r 
transformations. 

For three flavors, the flavor operators ta also denote the three antisymmetric Gell- 
Mann matrices, i.e., A = 2,5, 7. In this case the two-flavor condensation pattern, Eq. (|2]) 
is still possible, but now there are several other combinations which cannot be transformed 
into S22 via color or flavor rotations. If we have three degenerate light flavors, dense matter 



is expected to form a so-called color-flavor locked (CFL) state ]T5[, characterized by the 
situation 

S22 = S55 = s 77 7^ and saa' = if A^A'. (3) 
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In this state color SU(3) as well as the chiral SU(3)l x SU(3)r and the [/(l)-symmetry 
related to baryon-number conservation are broken down to a common SU(3) co i or+ y sub- 
group where color and flavor rotations are locked. 

Both situations discussed so far are idealizations of the real world, where the strange 
quark mass M s is neither infinite, such that strange quarks can be completely neglected, 
nor degenerate with the masses of the up- and down quarks. However, for sufficiently large 
quark chemical potentials /i 3> M s , the strange quark mass becomes almost negligible 
against fi and a phase very similar to Eq. (|3]) should form. In this case one expects that 
S55 and S77 are somewhat smaller than S22 but do not vanish. Moreover, if we assume 
an unbroken isospin symmetry, M u = Mj, which we will always do in this article, we 
still should have s 55 = s 77 . This corresponds to an SU{2) co i or+v subgroup of the original 
symmetry, where isospin rotations are locked to certain SU(2) rotations in color space. 
Therefore this phase is usually called CFL-phase as well. 

At low chemical potentials, the strange quark mass cannot be neglected against \i and 
a phase with s 55 = S77 7^ is no longer favored. On the other hand, below a certain value 
of n, one should finally reach the hadronic phase. Hence, if we start in the CFL-phase 
and keep lowering the chemical potential, there are (at least) two possible scenarios: (1) 

(2) 

Below a critical value jj c , S55 and S77 become zero, but S22 remains non-zero, i.e., we have 
basically a transition into the 2SC phase (with or without uncondensed strange quarks). 
Upon further decreasing fi we finally reach the hadronic phase at some critical chemical 
potential fi^ . (2) The CFL-phase is directly connected to the hadronic phase without a 
2SC-phase in between. This scenario has the interesting feature of the so-called quark- 
hadron continuity, which is related to the fact that the symmetries in both regimes are 
the same and the low-lying spectra in both phases can be mapped onto each other after 
the appropriate redefinition of the conserved U(l) charge in the CFL-phase [|lf|. 

It is obvious that the answer to the question which of the two scenarios is realized 
in nature depends on the strange quark mass. This has first been analyzed by Alford, 
Berges and Rajagopal [|T7| who have studied the color-flavor unlocking phase transition 



in a model calculation with different values of Ms- Assuming that the region below 
fi ~ 400 MeV belongs to the hadronic phase, these authors came to the conclusion that 
a 2SC-phase exists if M s > 250 MeV. Here M s is an effective "constituent" mass of the 
strange quark, which could be considerably larger than the current quark mass m s ~ 100 
to 150 MeV in the Lagrangian. In general, it is T- and /i-dependent. Moreover, it could 
depend on the presence of quark-antiquark and diquark condensates. In particular, it can 
be discontinous along a first-order phase transition line. This means, not only the phase 
structure depends on the effective quark mass, but also the quark mass depends on the 



phase. This interdependence has not been taken into account in Ref. flT7|, where M s has 



been kept at fixed values. The authors of Ref. Po] , who have studied the effect of a 
nonzero strange quarks mass within the framework of an instanton mediated interaction, 
also neglected these interdependencies. 

The T- and /i-dependence of effective quark masses in non-color superconducting 
phases has been studied extensively within NJL-type models nT9L EO, EHI], where the masses 



are closely related to the (T- and /z-dependent) quark-antiquark condensates 

4> u = (uu) = (dd) and <fi s = (ss) . (4) 
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Since color-superconductivity is also very often studied within NJL-type models, it is the 
obvious next step to describe both, diquark condensates and quark-antiquark condensates, 
simultaneously within the same model. For two flavors this has been done by several 
authors with different degrees of sophistication |23|, ^4]. For three flavors, a first 
attempt in this direction was made in Ref. [p5|] . However, as we will discuss in the next 
section, the thermodynamic potential presented in that article does not have the correct 
behavior in the limiting case of a CFL-phase with exact 5'f/(3)-symmetry and does also 
not reduce correctly to the 2SC thermodynamic potential, when S55 and S77 are switched 
off.Q In the present article we therefore revisit the problem. In Sect. |2| we present our 
model and derive the thermodynamic potential in mean-field approximation. Our result 
has the correct behavior for the 2SC case and for the 5'L r (3)-symmetric CFL phase. In 
general, the dispersion laws of the corresponding quasiparticle states are more complicated 
than the result presented in Ref. |25|. In Sect. [|we present numerical results for zero and 
no n- zero temperatures. Finally, conclusions can be found in Sect. 0L 



2 Formalism 

We consider the effective Lagrangian 

£eff = Tpiifi-m)^ + C qq + C qq , (5) 

where tp denotes a quark field with three flavors and three colors. The mass matrix rh has 
the form rh = diag(m u ,md,m s ) in flavor space. Throughout this paper we will assume 
isospin symmetry, m u = rrid- 

Since we wish to study the interplay between the color-superconducting diquark con- 
densates saa' and the quark-antiquark condensates (p u and <p s , we need an interaction 
which allows for condensation in these channels. To this end we consider an NJL-type 
interaction with a quark-antiquark part 
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£qq = Gj^ [OKVO 2 + (M5TaV) 2 ] (6) 



a=0 



and a quark-quark part 

C qq = H J2 (i>il5T A \ A ,Ctlj T )(^ T Ct l5 T A \ A ,ij). (7) 

A=2,5,7 A' =2,5,7 

These effective interactions might arise from some underlying more microscopic theory 
and are understood to be used at mean-field level in Hartree approximation. 

In Eq. (|6]), r a ,a = 1,...,8, denote Gell-Mann matrices acting in flavor space, while 

r ° = yl ■"■/ ^ s P ro P° r ti° na l to the unit matrix. Hence C qq is a U(3)l x Z7 (3) ^-symmetric 
4-point interaction. In NJL-model studies of the meson spectrum [p6|, |27| or of properties 



of quark matter at finite densities or temperatures [28, 31], usually a 't Hooft-type 



* In the meantime a substantially revised version of Ref. |^5| appeared, see footnote on page 0. 
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6-point interaction is added which breaks the the £7a(1) symmetry. It is straight forward 
to take it into account such a term, but for simplicity we neglect it here. As we will 
discuss below, this might have qualitative consequences for the results. 

In order to calculate the mean-field thermodynamic potential at temperature T and 
quark chemical potential [i, we linearize the effective Lagrangian in the presence of the 
expectation values (fi u = (fid, 4>s, $22, and s 55 = s 77 . In this context it is convenient to 
introduce the effective quark masses 

M u = m u - AG(fi u and M s = m s - AG(fi s (8) 

and the diquark gaps 

A 2 = -2Hs 2 2 and A 5 = -2Hs 55 . (9) 
Then, after formally doubling the degrees of freedom, 

:= ± (SL) , HO) 



y/2 \Clfi-{X, y 

the thermodynamic potential per volume can then be written as 

Q(T,fi) = -Tj^J-^ ^r In (^S~\iu n ^\ 

+ AG<fil + 2G<fi 2 s + H(\s 22 \ 2 + 2| S55 | 2 ). (11) 

The inverse propagator of the g-fields at 4-momentum p is given by 

( + A 275 r 2 A 2 + A 5 7 5 (r 5 A 5 + r 7 A 7 ) \ 

KP) \-A^ 5 t 2 \ 2 -A* 515 (t 5 \ 5 + t 7 A 7 ) 1 ' 

with M = diag(M u , M u , M s ). In Eq. (0), S _1 has to be evaluated at p — (iu) n ,p) where 
u n = (2n — 1)ttT are fermionic Matsubara frequencies. 

Taking into account the Dirac structure, color, flavor, and the charge conjugate com- 
ponents, the g-fields are 72 dimensional objects, and the trace in Eq. fllTD has to be 
evaluated in this 72 dimensional space. A tedious but straight-forward calculation yields: 

^Tr In ^-S ,-1 (iw n ,^) 

= 31n(ij(x+ M x^ + 2|A 2 |V« + |A 2 | 4 )) 

+ 21n(ij(x+x" +2|A 5 |V S +|A 5 | 4 )) + 2 In (^(z+z" + 2| A 5 \ 2 y us + | A 5 | 4 )) 

"I" ^ n (^jts ( x uu x uu x ss x ss ~l~ 2| A 2 | x ss x ss y uu + 4| A5I {oc uu x ss + x ss x uu )y us 
+ \A 2 \ 4 x+x~ + 4|A 2 | 2 |A 5 | 2 (x+x- + x+ s x~) 
+ M^5\ 4 (x^ u x~ + x+x~ u + Ay 2 us ) 

+ 8| A 2 | 2 |A 5 |V + 32|A 5 | 6 y us + 16|A 5 | 8 ) ) , (13) 
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where we introduced the abbreviations 

„± I, . _L „-..\2 i ^2 



X 



ff 



(u> n ±ifi) + p + MfMfi and ?//// 



ul + fi 2 +p 2 + M f M r . (14) 



With these definitions one finds for the argument of the first logarithm on the r.h.s. of 
Eq. (m 



+ 2\A 2 \ 2 y uu + |A 2 



(o£ + £- 2 )(^ + 3f) 



with 



(15) 



(16) 



These are exactly the dispersion laws of the paired quarks in a two-flavor color supercon- 
ductor and the corresponding Matsubara sums are readily turned out using the standard 
relation 



-|A fe |/TN 



(17) 



The other terms in Eq. (|13|) are in general more complicated. There are, however, two 
simplifying limits. The first one corresponds to a two-flavor color superconductor, together 
with unpaired strange quarks. In this case A5 vanishes and Eq. ([13]) becomes 



-T: In (J-S l (iu n ,p) 
,ul + E- 2 



-- 4 
+ 3 



ln( 
ln( 



Ji2 
J^2 



A 5 =0 

1 j 2 4. W+ 2 



2^2 

uj 2 +et 2 . 



+ 2 



ln(^ 



wl +e, 2 



2*2 



, , 2 4- P+ 2 
' UJ n + £ u - 

J^2 



J^2 



18) 



with £-, 



'f ~ yP 2 5 ^~ ^' Here we recover the fact, that only four of the six light quarks 
(two colors) participate in the 2SC condensate, while the two remaining ones (and of 
course all strange quarks) fulfill the dispersion laws of free particles with effective masses 
M f . 

We can also reproduce the structure of the dispersion laws of the idealized 3-flavor 
symmetric CFL-state. To this end we evaluate Eq. flilf ) for M u = M s and A2 = A5. One 
finds 



^Tr In (^-S \iuj n ,p) 



M U =M S , A 2 =A 6 



ln( 



ul + E, 



act 



J>2 



+ ln( 



w + E, 



act 



J>2 



+ 



.2 1 771- " 
f2~^- 



+ E si{ 



sing 



J^2 



(19) 



with Eg Ct = E± and Et 



sing 



(a/p 2 + M 2 ± /i) 2 + |2A 2 | 2 . We see that in this limit 
all quarks participate in a diquark condensate, forming an octet with diquark gap A 2 
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(= A5) and a singlet with diquark gap 2A2. In fact, it has been shown that for three 
degenerate flavors, the diquark gaps form a singlet with A sing = 2A§ + 4A 6 and an octet 
with A oct = A3 — Ag, where 3 and 6 refer to pairing in the color-antitriplet and color- 
sextet channel, respectively Jl5|, |32], |33|. Since we have neglected the small color-sextet 



contribution in our ansatz, we should find A sing = 2A orf , in agreement with the above 
results. 

Both limiting cases, Eqs. (|18|) and (fL9|), are not correctly reproduced, if one starts 
from the thermodynamic potential given in Ref. P^| . For instance, there are always 
quarks which fulfill free dispersion laws (no diquark gaps), which should not be the case 
in the CFL-phase, see abovej] 

The Matsubara sums over Eqs. (|i8|) and ([19|) can again be turned out with the help of 
Eq. (|TT|). In general, i.e., for Eq. ( pT3[ ) with arbitray values of the condensates, this cannot 
be done so easily. If one combines the second with the third logarithm on the r.h.s., the 
argument becomes a polynomial of fourth order in The same is true for the argument 
of the fourth logarithm. The corresponding dispersion laws are related to the zeros of 
these polynomials. Although, in principle, the zeros of a polynomial of fourth order can 
be determined analytically, the results are very nasty expressions which are difficult to 
handle. Therfore, in practice one has to determine the dispersion laws numerically. After 
that, one can again employ Eq. (|T7|) to calculate the Matsubara sum. Alternatively, 
one can turn out the Matsubara sum numerically without previous determination of the 
dispersion laws. This is the method we have used. 

So far, the thermodynamic potential depends on our choice of the condensates <fi u , 
4>s, S22, and S55. On the other hand, in a thermodynamically consistent treatment the 
condensates should follow from the thermodynamic potential by taking the appropriate 
derivatives. The self-consistent solutions are given by the stationary points of the poten- 
tial, 

5tt 5tt 5tt 5tt 5tt 5tt 

= 0. (20) 



5<p u 5<p s Ss 2 2 8si, 2 Ss 55 Ss* 55 
In this way one finds that the quark-antiquark condensates are given by 

d 3 p d 



<t>u = — T / - — ^77 — — Tr In ( — S 1 (iuJ n , p) 

YU 2 „ i O) dM u 2 VT V n,y) , 

<Ps = -TY / 7^ J- -Tr In (-S-\tuj n ,p)). (21) 
^ s ^ J 2tt 3 dM s 2 VT V n,FJ ) K ' 

The explicit evaluation of the derivatives is trivial, but leads to rather lengthy and not 
very illuminating expressions, which we do not want to present. Note, however, that M u 



t After submission of our final manuscript, Ref. |E5j was replaced by a substancially revised version 
(F. Gastineau et al, \iep-ph/010128!\ v3). The quark dispersion laws presented there are now consistent 
with the limits Eqs. (|18|) and (|19[). The authors make the simplifying assumption that the particle part 
and the antiparticle part of the Hamiltonian separate. This enables them to derive relatively simple 
analytical expressions for the quasiparticle energies. Although in general these solutions are not exact 
(corresponding to the zeros of the polynomials in our Eq. (|l3|)), the numerical results obtained in this 
way look very similar to ours, thereby justifying the above approximation. 
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and M s , which determine the right hand sides, depend again on <p u and <p s (see Eq. (§)) 
and therefore the equations have to be solved self-consistently. 

Similarly one finds that the diquark condensates are governed by the gap equations 

A 5 = 2i,A 5 r£ / (gi ^ 5* (is-'^D) ■ (22) 

Obviously, both equations have trivial solutions = which are independent of each 
other. For the nontrivial solutions Eqs. ([21] ) and (|22|) are coupled and have to be solved 
simultaneously. If there is more than one solution, the stable solution is the one which 
corresponds to the lowest value of Q. 



3 Numerical results 

For the numerical evaluation of the equations derived in the previous section we first have 
to specify the interaction. As a typical example we consider a four-fermion interaction 
with the quantum numbers of a single-gluon exchange, 

8 

C int = -g J2(^ X ^) 2 - ( 23 ) 

a=l 



This interaction was also the starting point of the model calculations in Refs. |L5[ and 



T7|. Performing Fierz transformations we obtain for the effective coupling constants 



which enter Eqs. (§) and (0) 



G = -g, H = -g. (24) 



In addition, C int gives rise to effective couplings in other channels, which in principle 
have to be taken into account to be fully self-consistent in the sense that also all exchange 
(Fock-) terms are included. In fact, already for two flavors, one finds that a self-consistent 
treatment requires the simultaneous consideration of possible expectation values in six 



different channels, like color dependent quark condensates or densities ||24|| . For three 
flavors, there are even more. We already mentioned the induced color-sextet diquark 
condensates in the CFL-phase [|l^, |32], |33| . However, most of these induced condensates 
- as far as they have been studied so far - have turned out to be small. Therefore, since 
the structure of the equations is rather involved even without these additional channels, 
we neglect them in this paper. Also, as stated above, Eq. (^3J) should only be viewed as 
a "typical" interaction, which we have chosen for a qualitative illustration of the more 
general equations derived in the previous section. Certainly, there is no reason to exclude 



other terms (see also the comments in Ref. |L7| about this point). 

Taking the model as defined above, there are four parameters: the coupling constant g, 
the bare quark masses m u and m s , and a cut-off parameter, which is needed to regularize 
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Figure 1: Gap parameters at T = as functions of the quark chemical potential fi. Left: 
Constituent masses of up and down quarks (dashed), and of strange quarks (solid). Right: 
Diquark gaps A2 (dashed) and A5 (solid). 



the divergent momentum integrals in Eqs. (]II), (PH) and (|22|). For simplicity, we take 
a sharp 3-momentum cut-off A. We expect that the results will be qualitatively the 
same, if a smooth form factor is used. As customary (see e.g. the comments on this 



point in Ref. |22[]), the parameters are chosen such that "reasonable" vacuum properties 
are obtained. We take A = 600 MeV, gA 2 = 2.6, m u = 5 MeV and m s = 120 MeV. 
With these parameters we find the vacuum constituent quark masses M u = 362.5 MeV 
and M s = 557.4 MeV, which corresponds to the condensates <p u = (— 240.5MeV) 3 and 
(p s = (— 257.3MeV) 3 . These values are similar to those obtained in Ref. f29| , where the 
parameters have been fixed by fitting vacuum masses and decay constants of pseudoscalar 
mesons. 

We begin with the discussion of the results at zero temperature. The behavior of the 
four gap parameters as functions of the quark chemical potential fi is displayed in Fig. [1|. 
In the left panel we show the constituent quark masses M u and M s , in the right panel 
the diquark gaps A 2 and A 5 . One can clearly distinguish three phases. At low chemical 
potentials \i < jMi — 342.8 MeV, the diquark gaps vanish and the constituent quark masses 
stay at their vacuum values. For the condensates this means s 22 = S55 = S77 = 0, while 
<p u and cf) s are large. Hence, in a very schematic sense, this phase can be identified with 
the "hadronic phase", although there are no hadrons in our model. In fact, in Fig. |I| the 
entire region below /i = /ii corresponds to the vacuum solution with zero density. 

At fi = fii a first-order phase transition takes place and the system becomes a two- 
flavor color superconductor: The diquark condensate S22 has now a non-vanishing ex- 
pectation value, related to a non- vanishing diquark gap A2, whereas A5 remains zero. 
Just above the phase transition we find A2 = 111.5 MeV. At the same time the mass of 
the up quark drops from the vacuum value to M u = 55.3 MeV. With increasing /z, M u 
decreases further, while A 2 increases until it reaches a maximum at fi ~ 475 MeV. Just 
below the next phase transition at fi — {12 — 508.5 MeV we find A 2 = 135.6 MeV and 
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M u = 8.6 MeV. 

In the 2SC phase the baryon number density is no longer zero and increases from 
Pb = 0.42 fm -3 (2.5 times nuclear matter density) above \x = p,\ to ps = 1.15 fm -3 
(6.7 times nuclear matter density) below p, = /z 2 . The density of strange quarks remains 
zero up to /i = /z 2 , but this could be an artifact of the interaction which contains no 
flavor mixing terms. As a consequence the strange quark mass M s remains at its vacuum 
value and hence strange quark states cannot be polulated in this regime. A flavor mixing 
interaction would lead to a reduced value of M s in the 2SC phase and thus to a lower 
threshold for a finite density of strange quarks. 

At p, = /j,2 the system undergoes a second phase transition, this time from the 2SC 
phase into the CFL phase, which is characterized by a non-vanishing diquark gap A 5 
(together with a non- vanishing A 2 ). The phase transition is again of first order: At the 
transition point A 5 jumps from zero to 109.7 MeV, while M s drops from 557.4 MeV to 
204.3 MeV. The non-strange quantities are also discontinous and change in the opposite 
direction: A 2 drops from 135.6 MeV to 111.2 MeV, and M u jumps from 8.6 MeV to 
17.5 MeV. The density jumps from p# = 1.15 fm~ 3 to 1.63 fm -3 , i.e., from 6.7 to 9.6 
times nuclear matter density. 

If we started from an exact SU(3) flavor symmetry, i.e., equal masses for all quark 
flavors, we would expect A 2 = A5 in the CFL phase. However, as we are far away from 
an exact symmetry, it is remarkable that the diquark gaps A 2 and A5 are so similar at the 
transition point. The reason for the phase transition is not that u and s quarks can no 
longer condense below that value. We still find a CFL solution down to p, = 445.0 MeV. 
However, for p, < /i 2 , there is a more favored solution with a much higher strange quark 
mass which cannot support a condensation of u and s quarks. In this regime the CFL 
solution corresponds to a metastable state. 

In Ref. JT7| it has been argued that the color-flavor-unlocking transition at zero tem- 
perature has to be first order because pairing between light and strange quarks can only 
occur if the gap is of the same order as the mismatch between the Fermi surfaces. Hence 
the value of the gap must be discontineous at the phase transition. The existence of pair- 
ing between light and strange quarks does of course not automatically mean that the CFL 
solution corresponds to an absolutely stable solution. This has already been emphasized 
in Refs. [[34], [35]. The quantitative criteria derived in these references for the existence of 
a CFL solution as a global minimum of the thermodynamic potential cannot be applied 
to our case because the possibility of a density dependent strange quark mass has not 
been taken into account in the derivation. However, at least qualitatively, our results 
support the arguments put forward in these references: There is indeed a non-vanishing 
lowest possible value of A5 for metastable CFL solutions and the minimal value of A5 in 
an absolutely stable CFL phase is even larger. 

At large p, we become more and more sensitive to the cut-off and we approach the 
limits of the model. This is also the reason why at large chemical potentials both diquark 
condensates decrease with p,. Although it is sometimes argued that the cut-off simulates 
asymptotic freedom, i.e., the weakening of the strong coupling constant at large (Fermi-) 
momenta in some crude way, the results should certainly not been taken too serious in 
this regime. In particular the fact that we find A5 > A 2 above \x ~ 520 MeV should be 
seen as a cut-off artifact. 
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300 400 500 
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Figure 2: Phase diagram in the — T plane. The solid lines indicate phase boundaries with a 
first-order phase transition, the dashed lines correspond to second-order phase transitions. The 
different phases can be distuingished by different values for the various condensates. Within 
the figure we have indicated only those condensates which are significantly different from zero. 
(Note that because of isospin symmetry, 4>d = <p u and sn = S55.) 



We now extend our analysis to non-vanishing temperatures. The resulting phase 
diagram in the fi — T plane is shown in Fig. 0. Let us first concentrate on the quark- 
antiquark condensates (f) u and <p s . We could then distinguish between three different 
regimes: At low temperatures and low chemical potentials both condensates, <p u and 4> s 
are relatively large and we can schematically identify this phase with the hadronic phase, 
as before. When we increase the temperature or the chemical potential, eventually a 
second phase is reached, where <p u is almost zero. S is still large in this phase and it 
becomes small only in a third phase at even higher temperatures or chemical potentials. 
At low temperatures, e.g., along the T=0-axis, these three phases are well separated 
by first-order phase transitions, whereas at high temperatures we only find smooth cross- 
overs. Consequently, both first-order phase transition lines end in a second-order endpoint, 
which in our model is located at T ~ 66 MeV and fi ~ 322 MeV for the phase transition 
related to <ft u and at T ~ 66 MeV and fi ~ 505 MeV for the phase transition related to 
<fi s . For the former we should keep in mind, however, that our model does neither contain 
hadronic nor gluonic degrees of freedom, which are certainly important in this area of the 
phase diagram. 

We now turn to the diquark condensates. As discussed for the zero temperature case, 
at low temperatures, the three phases discussed in the previous paragraph coincide with 
the "hadronic phase", the 2SC phase and the CFL phase, i.e., we have s 2 2 = s 55 = 
at low chemical potentials, S55 = 0, but S22 7^ at intermediate chemical potentials, 
and S22 7^ and S55 7^ at high chemical potentials. At high temperatures all diquark 
condensates vanish. This leads to two additional phase transition lines. The first one 
corresponds to the transition from the 2SC phase to the "quark-gluon" plasma, i.e., to 
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Figure 3: Gap parameters at fi = 520 MeV as functions of the temperature. Left: Constituent 
masses of up and down quarks (dashed), and of strange quarks (solid). Right: Diquark gaps A2 
(dashed) and A5 (solid). 



the phase with vanishing diquark condensates and (almost) vanishing <p u . In agreement 
with various BCS-type calculations, this phase transition is of second order. In addition 
we obtain almost perfect agreement with a universal relation which can be derived (see 
e.g. IllJ) between the value of the gap parameter (in our case this corresponds to A 2 ) 



at zero temperature and the critical temperature T c for that transition: 



T c « 0.57A 2 (T = 0) . 



(25) 



For instance, at fj, = 343 MeV we find A 2 (T = 0) = 111.5 MeV, and T c = 62.0 MeV, 
whereas from the above relation we would expect T c = 63.5 MeV. Since Eq. (25) is an 
approximate relation, based on the assumption that pairing takes place only in the vicinity 
of the Fermi surface, the agreement is quite good. 

Starting from the CFL phase and increasing the temperature, one first observes a 
"melting" of the diquark gap A 5 , before at somewhat higher temperatures A 2 also van- 
ishes. The intermediate 2SC phase "above" the CFL-phase is distinguished from the 
regime "left" to the CFL-phase by a much smaller value of S and hence of the strange 
quark mass. Based on the assumption that A5 is always smaller than A 2 (which is not 
really true at T = 0, as we have seen), the earlier disappearance of A5 and hence the exis- 
tence of a 2SC phase "above" the CFL phase was already anticipated in Ref. ||. Applying 
similar arguments as at zero temperature (cf. Refs. |T7|, [34], Q) it was also predicted in 
that reference, that the corresponding color-flavor unlocking transition should again be of 
first order. This was corrobated by a second argument, claiming that the phase transition 
corresponds to a finite temperature chiral restoration phase transition in a three-flavor 
theory. In this case the universality arguments of Ref. [3D] should apply, stating that the 
phase transition should be of first order. 

Indeed, following the phase transition line from the left, we find that the transition 
continues to be first order. However, above a tricritical point at fi ~ 518 MeV and 
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T ~ 59 MeV the phase transition becomes second order. This is illustrated in Fig. [^, 
where the constituent masses (left panel) and the diquark gaps (right panel) are displayed 
as functions of the temperature for fixed p = 520 MeV. We find a second-order phase 
transition at T = 59.8 MeV. Of course we cannot exclude that the second-order phase 
transition is just an artifact of our simplified model or the present approximation. On the 
other hand, the arguments of Refs. fT7| , [34], |35| in favor of a first-order phase transition 
are much less stringent at finite temperature, where even for vanishing condensates the 
Fermi surfaces are smeared due to thermal effects. The applicability of the universality 
argument is also questionable in the present situation because the 2SC phase is not a 
three-flavor chirally restored phase, but only SU(2) x SU (2) symmetric. Note that even 
SU(3)v is broken explicitly by unequal current quark masses and spontaneously, e.g., 
by the flavor-antitriplet diquark condensates. Therefore a rigorous prediction of the true 
order of the phase transition is rather difficult. 

Finally, we should repeat that our results at p > 500 MeV are likely to be sensitive 
to the cut-off. Therfore it remains to be studied in a more systematic way, which of our 
findings persist if the model parameters (including the cut-off) are varied over a wider 
range. Work in this direction is in progress. 



4 Conclusions 

Recent investigations revealed a rich phase structure of strongly interacting matter. Be- 
sides the various known phases of hadronic matter and the quark-gluon plasma phase at 
high temperatures several "superconducting" phases are likely to exist in cold dense mat- 
ter. Calculations within effective models showed that in this context "dense" could mean 
only a few times nuclear matter density such that those superconducting phases could be 
relevant for the interior of neutron stars or even for relativistic heavy ion collisions ||. In 
this paper we focused on the study of this "cold dense" region within the framework of 
an NJL-type effective quark model. This model allows for condensation in various chan- 
nels, in particular chiral symmetry breaking due to quark- ant iquark condensation and 
diquark condensation can be described. Our main intention was to examine the effect 
of selfconsistently (via gg-condensation) determined realistic quark masses. In particular 
the mass of the strange quark plays an important role. For the two extreme cases, taking 
the strange quark mass equal to the light quark masses or considering an infinite strange 
quark mass, two different condensation patterns are found, respectively: The so-called 
color-flavor-locked phase Hl5| , where all quarks contribute to diquark condensation and 



the two-flavor color superconductor || [| , where only light quarks participate in the con- 
densation. We could show that both limiting cases are correctly described within our 
model. However, as in the regions of interest the strange quark mass is neither negligi- 
ble nor infinite the realistic situation must be somewhere in between. Indeed, for quark 
matter at zero temperature, we find that for moderate values of the chemical potential 
the two-flavor superconductor is energetically favored, corresponding to densities of about 
2.5 po U P to about 6.7 po- Only for even higher densities the CFL phase is found to be 
the energetically most favored one. The main reason is the large constituent mass of the 
strange quarks compared with that of up and down quarks which has been determined 
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selfconsistently within our investigations. In particular, the strange quark mass in the 
2SC phase is much larger than than in the CFL phase and the phase transition is mostly 
triggered by a discontinous change of the quark mass. 

The situation remains qualitatively unchanged up to temperatures of about 50 MeV. 
At higher temperatures the diquark condensates melt. For the diquark condensates which 
involve strange quarks this happens at somewhat lower temperatures than for the non- 
strange condensate. Hence, with increasing temperature, the system passes from the 
CFL phase through an intermediate 2SC phase to the quark-gluon plasma. In contrast 
to the situation at zero temperature, we found that the strange quark mass in this "high- 
temperature" 2SC phase is comparable to the strange quark mass in the CFL phase. 
Furthermore we found a tricritical point on the color-flavor unlocking line above which 
the phase transition is second order. 

Of course, our model is still rather schematic and can only give some qualitative hints 
on the interdependencies between the effective quark masses and the various supercon- 
ducting phases. Unfortunately, as our knowledge about the effective quark interaction 
to be used at moderate chemical potentials is limited, more quantitative predictions are 
difficult. In this situation a systematic examination of the sensitivity of the results on 
the choice of the interaction might be helpful to sort out some common features. Perhaps 
the most severe simplification in our model was that we did not include a 't Hooft-type 
interaction, which mixes different flavors and thus affects the strange quark mass even in 
regimes, where no strange quark states are populated. It is therefore certainly worthwhile 
to investigate how such a term influences the phase structure. In a more systematic anal- 
ysis one should also study the effect of other condensates, which have been neglected so 
far. 

For simplicity, we restricted our studies to a common baryon chemical potential for all 
flavors. However, for many applications, like quark cores of neutron stars or the possible 
existence of absolutely stable strange quark matter, it would be interesting to consider 
/3-equillibrated quark matter, where the chemical potential for up quarks is in general 
different from that of down and strange quarks. This could alter the position of the 
various Fermi surfaces and hence the entire phase structure. In an earlier work without 
diquark condensation we have shown that a selfconsistent treatment of the strange quark 



constituent mass is a very important effect in this context pif. On the other hand, in 



more recent studies which include diquark condensation (e.g., about the existence of the 
color-flavor locked phase in neutron stars |36j or strangelets |37| ) the interdependencies 
between strange quark mass and phase structure have so far been neglected. In view of 
our present results, work in this direction seems to be worth being further pursued. 
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